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Human biomonitoring studies in the Canadian Arctic liave measured a wide range of metals and persistent 
organic pollutants in Aboriginal and non-Aboriginal mothers during two time periods in the Northwest 
Territories and Nunavut. This analysis provides preliminary estimates on sample sizes and sampling 
frequencies required to measure significant changes in maternal blood concentrations for PCB 153 and total 
mercury. For example, sample sizes of 35^0 mothers permit the detection of a 40% decrease in these chemical 
concentrations between two groups (e.g. communities or regions). Improvements in method sensitivity can be 
achieved by on-going sampling over multiple time periods (e.g. 4 or 5) in these regions, or increasing sample 
sizes. 
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Persistent organic pollutants (POPs) and mercury 
biomagnify in northern traditional foods and have 
been linked to human health effects, with most 
concern placed on the immune, nervous, hormonal or 
cardiovascular systems of the foetus. Accordingly, human 
biomonitoring studies conducted in the Canadian Arctic 
have frequently examined expectant women to obtain 
insight on potential risks to the developing child. 

Northern studies indicate that concentrations of 
many historic POPs have declined in Arctic biota ( 1 ) and 
Canadian Arctic peoples (2) over the past few decades. The 
fate of mercury in the Arctic is less clear, as both increasing 
and decreasing concentrations have been observed in biota 
(3), while human concentrations in Canada have shown a 
general decline (2). In view of this ambiguity, it becomes 
desirable to characterize temporal and spatial trends 
using clear quantitative measures. A statistically robust 
regression-based analysis method has been developed to 
describe temporal trends of POPs (1) and mercury (3) in 
Arctic biota. Although this approach has been undertaken 
for the general Inuit population in Greenland (4), to our 
knowledge this kind of statistical rigor has not been used 
to identify chemical trends or define sampling criteria for 
human biomonitoring in the Canadian Arctic. To inform 
future study design, we performed a preliminary statistical 



analysis of two previous northern biomonitoring studies 
to estimate sample sizes and sampling frequencies re- 
quired to report significant changes in blood concentra- 
tions for PCB 153 and mercury in mothers from the 
Northwest Territories and Nunavut. 

Methods 
Study population 

The Northern Contaminants Program (NCP) of Abori- 
ginal Affairs and Northern Development Canada has 
coordinated maternal biomonitoring studies in the Inuvik 
Region of the Northwest Territories and the Baffin Region 
of Nunavut on two occasions (1997-1999 and 2005-2007) 
(2). All expectant Aboriginal and non- Aboriginal mothers 
who volunteered for the NCP studies were sampled due to 
small population sizes, summarized in Table I. Signed 
informed consent was obtained from each participant. 
The study protocols were reviewed and approved by the 
research ethics boards for each of the participating centres 
and health authorities, as appropriate. 

Statistical analysis 

The polychlorinated biphenyl congener PCB 153 and total 
mercury were chosen for subsequent analysis because they 
were nearly 100% detected (99.04%) and they had amongst 
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Table I. Demographic variables for mothers from the Inuvik Region of the Northwest Territories, and the Baffin Region of Nunavut 



Inuvik 



Baffin 



Baseline study: 1 998-1 999 Follow-up study: 2005-2006 Baseline study: 1997 Follow-up study: 2005-2007 



Age (years) 

Mean (range) 
Ethnicity 

Inuit 

Dene/Metis 
Non-Aboriginal 
Parity (No. of children) 
1 
2 
3 

4 + 



25.6 (15-45) 

32.6% (n=31) 
44.2% (n=42) 
23.2% (n=22) 

40.4% (n=38) 
24.5% (n=23) 
19.1% (n = 18) 
16.0% (n = 15) 



25.3 (16-39) 

68.4% (n=54) 
24.1% (n = 19) 
7.6% (n=6) 

36.7% (n=29) 
22.8% (n = 18) 
12.7% (n = 10) 
27.8% (n=22) 



25.3 (15-39) 

88.6% (n=31) 

N/A 
11.4% (n=4) 

34.3% (n = 12) 
14.3% (n=5) 
17.1% (n=6) 
34.3% (n = 12) 



24.1 (15-39) 

100% (n = 101) 
N/A 
N/A 

21.8% (n=22) 
23.8% (n=24) 
24.8% (n=25) 
29.7% (n=30) 



the highest variabihty of all chemicals examined. Here, 
PCB 153 is expressed on a wet weight basis in plasma 
(|ig/L); total mercury is presented in whole blood (|ig/L). 
Probability plots and the Anderson-Darling test demon- 
strated that both chemicals were lognormally distributed. 
Statistical inferences were performed on log-transformed 
data; hence geometric means were assessed on the original 
scale for a multiplicative effect by determining the per- 
centage increase or decrease in chemical concentration 
that was detectable and significant. 

We first performed a temporal analysis where a 
population is re-sampled in subsequent time periods, 
making the assumption that each time point selects an 
independent sample of expectant mothers. While inde- 
pendence is typically not assumed for time series or 
longitudinal studies, the type of temporal dependence is 
difficult to determine from only two time periods. Further, 
the sampling design is such that each time period may 
necessitate a different sample of pregnant women, since 
the overall population is small and pregnancy is a 
temporary condition that is difficult to predict. We based 
these calculations on powers from the analysis of variance 
(ANOVA) hypothesis test. We note that these calculations 
were approximate since we assumed the estimate of the 
variability (mean square error or MSB) was constant for 
two or more time periods. If the MSB were to change, the 
calculated sample sizes would be different. 

We also examined a situation where 2 groups (e.g. 
region, ethnicity, or community) of equal size are sampled 
during the same time period. From the lognormal dis- 
tributions, we were able to declare the geometric mean 
of the first group as being significantly different from 
the second if it is greater than the corresponding upper 
confidence limit, or less than the lower confidence limit. 
This method of calculating confidence intervals on the log 
scale ensures the lower bound is always positive. As a 
result, however, they are no longer symmetric about the 



geometric mean. These group sample size calculations are 
based on the observed variation of a single pooled dataset 
for each chemical, expressed as the log of the standard 
deviation (logsd). 

Sample size calculations were based on 80'/o power 
(P =0.20) and using a significance level (a) of 0.05, at a 
95% confidence interval on the log scale. All statistical 
analyses were performed using the software package SAS 
Bnterprise Guide 4.2 (Statistical Analysis System). 

Results and discussion 

Table II presents our temporal analysis in terms of 
measurable percent decreases in geometric mean concen- 
trations for PCB 153 and total mercury, in anticipation of 
future decreases in human chemical concentrations. We 
observed that fewer samples are required as the number 
of time periods increases. For instance, in order to detect 
an overall 35% decrease in PCB 153 concentration, we 
would need to sample 40 mothers over each of 2 time 
periods, 33 mothers over each of 3 time periods, and so 
on. The number of samples required to measure the same 
percent change in total mercury is slightly higher. 

Our group analysis in Table III describes the percent 
change that can be measured between 2 populations 
during a single time period. For example, a sample size of 
40 participants (for each group) allows a 39.1% decrease 
or a 64.3% increase in the PCB 153 concentration to be 
detected. In other words, if the geometric mean concen- 
tration of the first group is 1 |ig/L, then the second group 
is significantly different if its geometric mean is smaller 
than 0.609 |.ig/L or larger than 1.643 p-g/L. As before, the 
estimates for total mercury are slightly larger. Chemicals 
with less variability will have a higher degree of precision 
using the sample sizes presented here. 

These results demonstrate that as study population 
sizes are increased, smaller percent changes in chemical 
concentrations can be detected. However, this approach 
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Table II. Percent decrease in geometric mean concentrations for PCB 153 and total mercury that can be detected for different sample 
sizes as the number of time periods increases (MSB =0.8757 for PCB 153; MSE = 1.0422 for total mercury) 



Log % Time Sample size Time Sample size Time Sample size Time Sample 

difference Decrease periods per period periods per period periods per period periods size per period 



PCB 1 53 










0.10 


10.5 


2 


345 


3 


0.15 


16.2 


2 


154 


3 


0.20 


22.1 


2 


87 


3 


0.30 


35.0 


2 


40 


3 


0.40 


49.2 


2 


23 


3 


0.50 


64.9 


2 


15 


3 


0.60 


82.2 


2 


11 


3 


0.70 


101 


2 


9 


3 


0.80 


123 


2 


7 


3 


0.90 


146 


2 


6 


3 


1.00 


172 


2 


5 


3 


lULdl lllclOUiy 










0.10 


10.5 


2 


410 


3 


0.15 


16.2 


2 


183 


3 


0.20 


22.1 


2 


104 


3 


0.30 


35.0 


2 


47 


3 


0.40 


49.2 


2 


27 


3 


0.50 


64.9 


2 


18 


3 


0.60 


82.2 


2 


13 


3 


0.70 


101 


2 


10 


3 


0.80 


123 


2 


8 


3 


0.90 


146 


2 


7 


3 


1.00 


172 


2 


6 


3 



283 


4 


240 


5 


210 


127 


4 


108 


5 


94 


72 


4 


61 


5 


54 


33 


4 


28 


5 


25 


19 


4 


16 


5 


15 


13 


4 


11 


5 


10 


g 


4 


8 


5 


7 


7 


4 


6 


5 


6 


6 


4 


5 


5 


5 


5 


4 


5 


5 


4 




A 


A 




A 


336 


4 


286 


5 


250 


150 


4 


128 


5 


112 


85 


4 


72 


5 


64 


39 


4 


33 


5 


29 


22 


4 


19 


5 


17 


15 


4 


13 


5 


11 


11 


4 


9 


5 


8 


8 


4 


7 


5 


7 


7 


4 


6 


5 


5 


6 


4 


5 


5 


5 


5 


4 


4 


5 


4 



for improving method sensitivity may present a challenge 
for future biomonitoring in the Canadian Arctic, parti- 
cularly during maternal studies, due to the inherent 
challenges of examining sparse populations in remote 
regions such as the Arctic (5). Improvements in method 
sensitivity may be better achieved by increasing the 



number of sampling time periods from 2 to 4 or 5 in 
order to achieve a precision of less than 20%. The 
northern data examined here indicates that it has been 
possible to confidently detect contaminant concentra- 
tion decreases in the range of 20 to 50'/o, depending on 
the chemical, based on prior sample sizes ranging from 



Table III Percent increase and decrease in geometric mean concentrations for PCB 153 and total mercury that can be detected when 
comparing 2 different groups of equal size (logsd = 1.132 for PCB 153; logsd = 1.203 for total mercury) 



PCB 153 



Total mercury 



Group sample size 


% Decrease 


% Increase 


% Decrease 


% Increase 


10 


62.9 


170 


65.2 


187 


15 


55.5 


125 


57.7 


137 


20 


50.4 


102 


52.6 


111 


30 


43.6 


77.4 


45.6 


83.8 


40 


39.1 


64.3 


41.0 


69.4 


50 


35.8 


55.9 


37.6 


60.2 


60 


33.3 


50.0 


35.0 


53.8 


70 


31.3 


45.5 


32.9 


49.0 


80 


29.6 


42.0 


31.1 


45.2 


90 


28.2 


39.2 


29.6 


42.1 


100 


26.9 


36.9 


28.4 


39.6 
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30 to 100 mothers (2). However, since these samples are 
not random population samples, care must be taken when 
making conclusions on sample size and temporal trends. 
We also acknowledge that with additional time points, 
other statistical approaches to analysing trends over time 
would involve regression and/or time series methods. 
Such methods are more computationally involved and 
typically require more than 5 time periods to adequately 
assess the changes in contaminant levels over time. 
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